$ cat input
3
1 2 3 1
5 6 2 2
4 7 5 3

$ ruby linear.rb input
Coefficient matrix and constant vector.
0) 1 2 3 1
1) 5 6 2 2
2) 4 7 5 3
Answer(x[i] & (A*x-b)[i]) follows
x[0]=-0.2666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666664e0  & 0.16e-99
x[1]=0.5333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333332e0  & 0.64e-99
x[2]=0.6666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666672e-1  & 0.4e-99

$ git switch master   
$ ruby linear.rb input
'bigdecimal/ludcmp' is deprecated and will be removed in a future release.
Coefficient matrix and constant vector.
0) 1 2 3 1
1) 5 6 2 2
2) 4 7 5 3
Answer(x[i] & (A*x-b)[i]) follows
x[0]=-0.2666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666664e0  & -0.2e-99
x[1]=0.5333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333333332e0  & 0.0
x[2]=0.6666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666666672e-1  & 0.0

{ans: [0.1218245836551854221294816349945599902708230426322897706646893441528370772984244758379821844214668379e1, 0.7182458365518542212948163499455999027082304263228977066468934415283707729842447583798218442146683795e0]}
{diff: [-0.1137205797e-98, -0.1372057968e-99]}